Displacement Echoes: Classical Decay and Quantum Freeze 
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Motivated by neutron scattering experiments, we investigate the decay of the fidelity with which 
a wave packet is reconstructed by a perfect time-reversal operation performed after a phase space 
displacement. In the semiclassical limit, we show that the decay rate is generically given by the 
Lyapunov exponent of the classical dynamics. For small displacements, we additionally show that, 
following a short-time Lyapunov decay, the decay freezes well above the ergodic value because of 
quantum effects. Our analytical results are corroborated by numerical simulations. 
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The fidelity with which a wavefunction is reconstructed 
after an imperfect time-reversal operation was originally 
introduced as a measure of reversibility in quantum me- 
chanics [lj]. Dubbed the Loschmidt Echo, it has re- 
ceived much attention in recent years in the context 
of decoherence and the quantum classical correspon- 
dence 0, [3, H, @, 0- For a generic perturbation of 
the Hamiltonian, four different decay regimes were found: 
the Gaussian perturbative regime, the Fermi Golden Rule 
(FGR) regime, the Lyapunov regime, and the regime of 
classically large perturbations. Of special interest is the 
Lyapunov regime where the purely quantum mechani- 
cal fidelity decays with the Lyapunov exponent of the 
classical dynamics. It suggests the existence of a uni- 
versal regime of environment-independent decoherence 
rate U&M- 

In this letter we analyze the decay of the Loschmidt 
echo under a new, non-generic perturbation, namely an 
instantaneous phase space displacement. Our investiga- 
tion is partly inspired by spectroscopies such as neutron 
scattering, Mossbauer 7-ray, and certain electronic tran- 
sitions in molecules and solids , [ll], EH ■ In these spec- 
troscopies phase space displacement (momentum boost 
or position shift) takes place with little or no change 
in the potential. Under this non-generic perturbation, 
we find that the decay rate of the average fidelity is al- 
ways set by the Lyapunov exponent. Moreover, for small 
displacements, the initial Lyapunov decay is followed at 
larger times by a quantum freeze of the fidelity at a 
displacement-dependent saturation value. Both scmiclas- 
sics and random matrix theory predict that the freeze 
persists up to infinitely large times. 

As our starting point, we recall that the differen- 
tial cross section for incoherent neutron scattering and 
Mossbauer emission/absorption can be calculated from 
the following correlation function 13|j 3 (from now on 
we set h = 1) 



Here, the brackets represent an ensemble average, fj are 
the position operators of the nuclei and H is the typical 
Hamiltonian of the target system. The ensemble average 
of the correlation function can be written [Io| : 
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where 



are coherent states with TV degrees of 



freedom, Q = Tr 



and $ (a) is a ther- 



mal weight, which tends to e P H ^( a > at high tem- 
peratures. The notation e lHpt = e~ jP ' r je jfft e lPr j 
suggests that we identify the kernel of the integral 
I(t) = He-'P ^e^VP ^e-^'la) with the kernel of 
a Loschmidt echo problem, we thus introduce the mo- 
mentum displacement echo 



M D (t) = \I(t)\ 2 = I (a\e 
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As introduction to our semiclassical calculation of the 
displacement echo Mo(t), we first discuss the validity 
of the diagonal approximation used in Ref. Q for the 
semiclassical approach to the Loschmidt echo and point 
out why this approximation is even better for the dis- 
placement echo. The diagonal approximation for the 
Loschmidt echo equates each classical trajectory 7 gen- 
erated by an unperturbed Hamiltonian H with a classi- 
cal trajectory generated by a perturbed Hamiltonian 
Hy = H + V. This procedure is not justified a priori in 
chaotic systems where one expects that an infinitesimally 
small perturbation generates trajectories diverging expo- 
nentially fast away from their unperturbed counterpart. 
It was however pointed out by Cerruti and Tomsovic [H , 
and later by Vanicek and Heller Q , that structural stabil- 
ity theorems justify this approximation. Roughly 
speaking one can show that, given a uniformly hyper- 
bolic Hamiltonian system H, and a generic perturbation 
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V, each classical trajectory -f' v generated by the (still hy- 
perbolic) perturbed Hamiltonian H + V is almost always 
arbitrarily close to one unperturbed trajectory 7. In gen- 
eral the two trajectories do not share common endpoints, 
however these endpoints are close enough that they are 
not resolved quantum-mechanically. This is illustrated in 
the left panel of Fig. 1. The semiclassical expression for 
the kernel of the Loschmidt echo involves a double sum 
over both the perturbed and the unperturbed classical 
trajectories, so that both and 7 are included. After a 
stationary phase condition, this double sum is reduced to 
a single sum where j v and 7 are equated. In other words, 
a semiclassical particle follows 7 in the forward direction, 
and 7y in the backward direction because this is the best 
way to minimize the action. The action difference is sim- 
ply given by the integral of the perturbation along the 
backward trajectory, and it is in general time-dependent. 




FIG. 1: Illustrative view of structural stability. Left panel: 
generic perturbation. 7 and 7' are two orbits of the unper- 
turbed Hamiltonian, yv is the orbit of the perturbed Hamilto- 
nian with the same initial condition as 7, while 7^ is the orbit 
of the perturbed Hamiltonian with the same initial condition 
as 7'. The endpoints of 7 and 7^ are separated by less than 
a quantum- mechanical resolution scale (shaded area). Right 
panel: phase space displacement. Labels are the same as in 
the left panel, with P replacing V as subscript for perturbed 
trajectories. Note that 7p and 7 lie on top of each other. 

In the case of a uniform phase-space displacement, the 
diagonal approximation becomes even better. This is 
so because any classical trajectory of the unperturbed 
Hamiltonian is also a trajectory of the perturbed Hamil- 
tonian, up to displacements at the trajectory's ends. This 
is illustrated in the right panel of Fig. 1. The fact that 
the action difference is here time- independent has the im- 
portant consequence that the FGR decay is replaced by a 
time-independent saturation term. The Lyapunov decay 
term is left almost unaffected, as it depends on the clas- 
sical measure of nearby trajectories with perturbed ini- 
tial conditions and does not depend on quantum action 
phases. We also note that for displacement echoes there 
is no Gaussian perturbative decay, since phase space dis- 
placements do not change the spectrum of the system 
aside from some possible but irrelevant global shift. 

For a quantitative approach to the problem, we scmi- 
classically evaluate Md (t) for the case of an initial Gaus- 



sian wavepacket, (r|a(r ,po)) = (7rg 2 ) 4 exp[zp • 
(r — r ) — (r — r ) 2 /2tr 2 ]. Following Ref. 17 1, we semi- 
classically propagate \a) with the help of the van Vleck 
propagator, linearly expanding around Tq: 



(r'\e- im \ a ) 
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^ y / C^exp[iS , 7 — mv^jl — a 1 (p 7 — p ) /2]. 



Here, the sum runs over all possible classical trajectories 
7 connecting r and r' in the time t, p 7 = —dS 1 /dr\ T=To 

S-j is the classical ac- 



is the initial momentum on 
tion accumulated on 7 

a 



7, 



v. 
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-0 2 Sj(r,r';i) /Ondr' 



is the Maslov index and 
. The kernel I(t) of 
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Mr>(t) involves a double sum over classical trajectories, 
7 and 7', and can be interpreted as the overlap between a 
wavepacket that is boosted and subsequently propagated 
with a wavepacket that is first propagated and subse- 
quently boosted Enforcing a stationary phase con- 
dition kills all but the contributions with the smallest 
actions. They correspond to 7 = 7' and one has 



I(t) 



v)V dr 'E elP(r ' +ro)c 7 
■t[c 1 
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x exp — — 1 (p 7 - po) 2 + (p 7 po - P) 2 



Taking the squared amplitude |/(£)| 2 one obtains the 
semiclassical expression for the displacement echo 
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We calculate (Mr>(t)}, the ensemble-averaged displace- 
ment echo over a set of initial Gaussian wavepackcts with 
varying center of mass ro. There are two qualitatively dif- 
ferent contributions to (Mj^(t)). The first contribution 
(M D (i)) c comes from pairs 7 m 7' of correlated trajecto- 
ries that remain within a distance < a of each other for 
the whole duration t of the experiment, while the second 
contribution (Md(£)) u arises from pairs of uncorrelated 
trajectories (7,7')- F° r t ne nrs t contribution, we write 
cxp[iP(r — r')] « 1, which is true in the semiclassical 
limit where a — > 0, and set 7 = 7'. One then has 

(M D (t)) c = ^drdrX(r-r') 

Q2 e -ff 2 [(p 7 -p ) 2 + (PT-P0-P) 2 ]\ (7) 
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where 5 a (r — r') restricts the integrals to |r — r'| < a. 
The calculation of (J7J) is straightforward. The integral 
over r' gives a factor a d . One then replaces one C 7 by its 
asymptotic value, C 7 oc exp[— Xt], and uses the second C 7 
to perform a change of integration variable J dr ^ 7 C 7 = 
/ dp. After a Gaussian integration, one finally gets the 
correlated contribution to (Md(£)) as 



(M D (t)) c = a exp[-(Pa) 2 /2] exp[-At]. 



(8) 



Here, a is a weakly time-dependent number of order 



For the uncorrclatcd part, an ergodicity assumption is 
justified at sufficiently large times, under which one gets 



(M D (t)) u = /(P) (M D (t)>«, 

/(P) = ^ / drdr ' exp[«P • (r - r') 



(M D (t)) u = 



(9a) 
(9b) 

(9c) 



x exp — - 



(p 7 - p ) 2 + (p 7 - p - P) 2 



with the system's volume il oc L d . It is straight- 
forwardly seen that (M-o(t)) n = cxp[— (Pcr) 2 /2], and 
/(P) = g(|P|L)/(|P|L) 2 , in term of an oscillatory func- 
tion g(\P\L) = 4sin 2 (|P|L/2) for d = 1 and g(\P\L) = 
4J 2 (|P|L) for d — 2. For d = 3, g is given by Bcssel and 
Struve functions. The uncorrclatcd contribution to the 
displacement echo reads 

(M D (i)) u = cxp[-(Pa) 2 /2]. 9 (|P|L)/(|P|L) 2 , (10) 

which, together with Eq. © gives the average displace- 
ment echo as 



(Af D (t)) = exp[-(Pa) 2 /2] 



.-At 



g(\P\L) 



(ii) 



In addition, as is the case for Loschmidt echoes, 
(Mu(t)) > N^ 1 where N is the size of Hilbert space. 

Eq. (fTTj) is our main result. It states that Md(£) is the 
sum of a time-dependent decaying term of classical ori- 
gin and a time-independent term of quantum origin. The 
latter can also be obtained within random matrix theory. 
The prefactor exp[— (Ptr) 2 /2] — ► 1 in the semiclassical 
limit and is thus of little importance. We see that gener- 
ically, Md (t) follows a classical exponential decay, possi- 
bly interrupted by a quantum freeze as long as the dis- 
placement is not too large, g(\P\L) j (|P|£) 2 > N^ 1 [3]. 

We note that in the semiclassical limit, Mu(t — ► 0) — ► 1, 
because of the saturation of a(t — ► 0) — ► 1 and the dis- 
appearance of uncorrelated contributions at short times. 
Most importantly, there is no displacement-dependent 
decay, i.e. no counterpart to the FGR decay nor the 



M D (t) 




FIG. 2: Main plot : Displacement echo Mn(t) for the kicked 
rotator model with N = 262144, and displacements P = 
m x 2n/N, m = 10, 20, 30. Averages have been performed 
over a set of 10000 different initial coherent states. The 
full lines correspond to kicking strengths K — 10.09, 50.09 
and 200.09 (from right to left) and the dashed lines (slightly 
shifted for clarity) give the predicted exponential decay given 
by the reduced Lyapunov exponent Ao = 1.1, 2.5, 3.7. The 
dotted line gives the saturation at A^ -1 . Inset : Displace- 
ment echo for A^ = 8192, K = 10.09, and displacements 
P = 2tt/N, 4tt/N, . . . Wtt/N. Data are obtained from 1000 
different initial coherent states. The dashed line gives the 
predicted exponential decay given by the reduced lyapunov 
exponent Ao = 1.1. The dotted line gives the minimal satu- 
ration value at A' -1 . 

perturbative Gaussian decay for Mrj(i), because phase- 
space displacements leave the spectrum unchanged, up 
to a possible homogeneous shift [l9l |. 

What does the "freeze" correspond to physically? It is 
the elastic component in any of the mentioned spectro- 
scopics: Mossbaucr, neutron, and molecular electronic, 
and was first identified by van Hove in connection with 
neutron scattering [l4j |. There is a finite probability, 
above the 1/N statistical limit, of not changing quan- 
tum state in spite of being "hit" ; this is the source for 
example of the recoilless peak in Mossbaucr spectroscopy. 

We now check our predictions numerically . We spe- 
cialize to the kicked rotator model with Hamiltonian 



P 

H = — + K cos x 8(t — n). 



(12) 



We focus on the regime K > 7, for which the dynamics is 
fully chaotic with Lyapunov exponent A = ln[if/2]. Wc 
quantize this Hamiltonian on a torus, which requires to 
consider discrete values pi = 2irl/N and xi = 2irl/N, I = 
1, ...N. In these units, one has L = N. The displacement 
echo of Eq. © is computed for discrete times t = n, as 

M D (n) = \(4>o\e- lP& (tf) n e iP& {U) n (13) 

with P = |P|. Here, we used the unitary Floquct 
time-evolution operator hi whose matrix elements in 
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FIG. 3: (color online) Main plot: Saturation value Md(oo) 
of the displacement echo as a function of the rescaled 
displacement NP/2n for the kicked rotator model with 
N = 256, 1024, 4096, 16384 (full lines, from top to bot- 
tom). Data are obtained from 1000 different initial co- 
herent states. The dotted lines give the saturation at 
N^ 1 . The red dashed line gives the theoretical prediction 

Md(oo) = Max(4exp[-( ( TP) 2 /2]sin 2 (PL/2) / /(PL) 2 ,Ar- 1 ) 

for TV = 16384. Inset: Quantum freeze of the displace- 
ment echo for kicking strength K — 10.09, N = 4096, and 
P £ [0,2n/N]. The dashed line gives the decay with the 
reduced Lyapunov exponent Aq = 1.1 (see text). 



x— representation are given by 



exp 



tt(Z-Z ') 
' N 



exp 



NK 2nV 

l cos 

2tt N 



The time-evolution of ip in Eq. (fT3")h is calculated by 
recursive calls to a fast-Fourier transform routine. 

Fig. 2 shows the behavior of the echo for displacements 
in the range P > 2tt/N for which (Mr>(t)) u < TV -1 
and thus plays no role. It is seen that the decay rate of 
the displacement echo strongly depends on the kicking 
strength K, but is largely independent of the displace- 
ment P. We quantitatively found that in that regime, 
M-£>(t) w expj— Apt], in term of the reduced Lyapunov 
exponent Ao [20(. The inset shows moreover, that low- 
ering the displacement to the regime P = m2n/N with 
m < 5 does not affect the decay rate of M D (i), i.e. there 
is no FGR decay for the displacement echo. 

We focus in Fig. 3 on smaller displacements P < 2ir/N. 
The behavior of (Md(£)) clearly satisfies Eq. (jlip . with 
a quantum freeze at a displacement-dependent value fol- 
lowing a decay with a slope given by the Lyapunov ex- 
ponent. We show in the main panel the P-dcpcndence 
of the value at which Mu(t) freezes. The data unam- 
biguously confirm the algebraically damped oscillations 
predicted in Eq. (fTTj) . 

In summary, we have presented a semiclassical calcula- 
tion of phase-space displacement echoes. We showed that 
they are generically given by the sum of a classical decay 
and a quantum freeze term (jlip . Because phase-space 
displacements do not generate time-dependent action dif- 



ferences, and because they vanish in first order pertur- 
bation theory, there is no other time-dependent decay, in 
contrast to Loschmidt echoes 

To conclude, we note that neutron scattering correla- 
tion functions are given by the average (I(t)) of the kernel 
of Mr>(t). Starting back from Eq. ([5]), one gets 



|F i ,(P,t)|^exp[-(P ( T) 2 /4] 



9 1/2 (\P\L) 
\P\L ■ 



(14) 



i.e. \Yjj \ is given by the quantum freeze term only. This 
is so, since the correlations between pairs of trajectories 
that are necessary for the existence of the Lyapunov term 
appear only once I(t) is squared. 
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